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Abstract 

We apply the dipole formalism that has been developed to describe low-x deep inelastic scattering 
to the case of ultra-high energy real photons with nucleon and nuclear targets. We hope that there 
will be future modeling applications in high-energy particle astrophysics. We modify the dipole 
model of McDermott, Frankfurt, Guzey, and Strikman (MFGS) by fixing the cross section at the 
maximum value allowed by the unitarity constraint whenever the dipole model would otherwise 
predict a unitarity violation. We observe that, under reasonable assumptions, a significant fraction 
of the real photon cross section results from dipole interactions where the QCD coupling constant 
is small, and that the MFGS model is consistent with the Froissart bound. The resulting model 
predicts a rise of the cross section of about a factor of 12 when the the photon energy is increased 
from 10 3 GeV to 10 12 GeV. We extend the analysis to the case of scattering off a 12 C target. We 
find that, due to the low thickness of the light nuclei, unitarity for the scattering off individual 
nucleons plays a larger role than for the scattering off the nucleus as a whole. At the same time 
the proximity to the black disk limit results in a substantial increase of the amount of nuclear 
shadowing. This, in turn, slows down the rate of increase of the total cross section with energy as 
compared to the proton case. As a result we find that the 12 C nuclear cross section rises by about 
a factor of 7 when the photon energy is increased from 10 3 GeV to 10 12 GeV. We also find that 
the fraction of the cross section due to production of charm reaches 30% for the highest considered 
energies with a 12 C target. 
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I. INTRODUCTION 


There is currently an interest in the types of showers induced by ultra-high energy (UHE) 
cosmic ray neutrinos |l| which will be relevant to the Anita, Auger and Icecube experiments. 
The resulting showers are presumed to be initiated by the Bremsstrahlung photons radiated 
from the electron produced in the initial reaction, 

v e -\- A —» e + X , (1) 


where A is a nucleus in the target medium, and X is a hadronic jet produced in the initial 


reaction. Due to the Landau-Pomeranchuk-Migdal (LPM) |2j effect, at UHE soft electro¬ 
magnetic radiation is suppressed and most of the energy of the electron is transfered directly 
to the photon. Furthermore, the cross section for e + e _ pair production drops at UHE and 
the hadronic interaction between the photon and the target nuclei may dominate while the 
electromagnetic interaction becomes negligible. The suppression due to the LPM effect be¬ 
comes stronger depending on the density of the target medium. For a general overview of 


the LPM effect and electromagnetic suppression various media, see Ref. |3[. Furthermore, 
the shapes of showers may depend upon whether they are dominantly electromagnetic or 


0 


hadronic |4||. Along these lines, it also important to determine what fraction of showers are 
due to charmed particles. This is important for IceCube and MACRO because an increase 
in the number of charmed mesons in the initial reaction will increase the number of high 
energy muons seen in experiments, but the contribution of high energy muons is one factor 
used to determine the composition of cosmic rays. In addition, charmed particles contribute 
to the flux of atmospheric neutrinos. Thus, experiments will need to take this contribution 
into account in searches for diffuse astrophysical neutrinos. For example, constraints on 


the UHE neutrino flux are sensitive (see, e.g., Ref. 5]) are sensitive to constraints on UHE 
photon cross sections. 

Another possible source of UHE cosmic photons is the decay of extremely masses exotic 


particles, topological defects, and Z-burst models [6[. Particles with masses as high as 10 26 eV 


Ref. 


0 


may explain the observation of super-GZK energy cosmic rays [7]. Indeed, the calculations in 


have shown, using both standard QCD and supersymmetric QCD, that a large part 
of the spectrum in the decay of super-massive particles consists of photons. A characteristic 
of these “top down” models is the existence of a large photon flux in cosmic rays. The ratio 
of protons to photons in the primary interaction for the production of showers can be used to 
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distinguish between various top-down scenarios. It is necessary, in order to address the issues 
above, to place upper limits on the growth of the real photon cross section. Upper limits 
have been placed on the the fraction of primary photons at 26% j^] by analysis at the Auger 
observatory, but those upper limits are sensitive to the photon-nucleus interaction. The 
main purpose of this paper is to combine the dipole picture (in the limit that the incident 
photon is real) with unitarity constraints and phenomenological expectations to estimate 
the actual growth of the real photon cross section with target nucleons/nuclei at UHE. We 
allow the cross section to grow with energy as fast as possible under the constraints of the 
S-channel unitarity, and thereby we at the very least place upper limits on the growth of 
photon-nucleus cross sections. It is worth emphasizing here that one cannot simply use a 
smooth extrapolation of the cross section to higher energies by assuming (as it is often done 
for the case of hadron-hadron scattering) a parameterization of the cross section inspired 
by the Froissart bound of the form, = a + Mn 2 (s/s 0 ). Asymptotically, the photon - 
hadron cross section can grow faster than the rate of growth supplied by the Froissart bound 
due to the the fact that the photon wave function is non-normalizable. In fact, as we will 
discuss, the rate of growth with energy of the photon-hadron cross section may be as fast as 
In 3 p, 43|. Because of the uncertainties involved in describing astronomically high energy 
photon behavior, the strategy that we take in this paper is to provide meaningful upper 
limits on the growth of the cross section based on the observation that the true hadronic 
cross section of the photon will be tamed in any realistic theory relative to what is predicted 
in a low order perturbative QCD description. Therefore, a low order perturbative QCD 
description that is only corrected for definite unitarity violations supplies a absolute upper 
limit. We purposefully avoid including a more precise description of higher order affects 
(using, e.g., saturation, color glass condesnsate, etc...) in order to avoid introducing any 
model dependent taming effects into the description of the high energy photon. Therefore, 
we most likely overestimate the cross section by a significant fraction, but we will see that 
the upper limits obtained are nevertheless useful since they provide stronger constraints 
than more complex models of the high energy behavior. 

Dipole models in strong interaction physics describe the interaction of a high energy 
virtual photon with a hadronic target by representing the photon by a distribution over 
hadronic Fock states of varying sizes in the light-cone formalism. This picture has become 
popular for studying deep inelastic scattering in the low-x limit. It is a dipole model because 
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the very small size configurations are usually modeled by a small quark/anti-quark pair, 
though larger size Fock states are not necessarily dipole configurations and should be thought 
of as more general hadronic states containing, for example, intrinsic gluon fields. The idea 
that hadrons should contain configurations interacting with different strength emerged first 
in the context of discussion of the inelastic diffraction Itj [ll, [l^ . This idea reemerged in 


connec 


ll M 


ion with the introduction of the two gluon exchange model for the strong interaction 


141 ]. In this model the interaction of a small hadronic system with a large hadronic 


system is obviously proportional to the transverse size squared of this system 


13]. In this 


context the impact parameter representation of the scattering amplitude first introduced for 


high energy processes by Cheng and Wu 3] in their studies of the high energy QED turned 
out to be useful. The relationship between the smallness of the interaction with nucleons, 
the small size of the incident configuration (color transparency), and Bjorken scaling was 


emphasized in 0 where it was demonstrated that small and large size configurations give 
comparable contributions to F 2 (x,Q 2 ) at small x and Q 2 ~ fewGeV 2 - the QCD aligned 
jet model. The use of the eikonal model with two gluon exchange and the Cheng and Wu 
representation for the impact parameter photon wave function (extended to finite Q 2 ) was 
used to build a model of nuclear shadowing jjj|. 


It was pointed out in [H| that the cross section of the dipole nucleon interaction within 
the leading log approximation is actually proportional to the gluon density at high virtuality 
and small x, resulting in a fast growth of the total cross section. The resulting dipole-nucleon 
cross sections were used to take into account finite Q 2 corrections to exclusive vector meson 


production in DIS 


0 


Splitting of the dipoles into systems of dipoles was studied within 
the BFKL approximation in j^oj]. The dipole model was applied to the description of HERA 


DIS data by Golec-Biernat and Wusthoff 121] who introduced a parameterization of the 
dipole cross section inspired by the eikonal model which ensured that the total cross section 
of the interaction of dipoles of any size would reach a finite limiting value at high energies. 
A more realistic model was introduced by McDermott, Frankfurt, Guzey and Strikman 


(MFGS) in Ref. 22], In |22f (the MFGS model) the cross section for small dipole sizes was 
constrained to satisfy the perturbative QCD expression for the dipole - nucleon interaction, 
while for large sizes a growth consistent with the pattern of the pion - nucleon interaction 
was imposed. An important advantage of the model is that one can adjust the behavior of 
large size and small size configurations independently. The MFGS model was extended to 
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an impact parameter analysis in Ref. 


earlier in Ref. 


23|, and other impact parameter analyses were done 


2J> 25]. 


At small Bjorken-x, the gluon distribution dominates over the quark distribution in the 
nucleon target, and the behavior of the small size dipole cross section is successfully predicted 
in leading-twist perturbative QCD (pQCD). However, at very small values of Bjorken-x, the 
leading-twist gluon distribution becomes unreasonably large and qualitatively new physics 
is expected to dominate. The usual assumption is that, whatever QCD mechanism is re¬ 
sponsible for extremely low-x behavior, scattering of a dipole of a given size at fixed impact 
parameter at sufficiently high energy occurs at or near the limit allowed by unitarity. Appli¬ 
cations of the dipole model have thus far been used mainly in searches for these qualitatively 
new regimes in QCD. Within the dipole picture, at fixed large Q 2 and sufficiently small x 
the distribution of sizes is sharply peaked around small size quark-antiquark pairs. For small 
photon virtualities, Q 2 , and for the energies E 1 < 10 3 GeV the distribution in sizes is dom¬ 
inated mostly by large sizes which lie far from the pQCD regime. However, for extremely 
high energy photons (relevant to cosmic rays), pQCD predicts a very rapid rise in the ba¬ 
sic cross section with an increase in dipole size. At asymptotically large photon energies, 
finite size configurations rapidly reach the unitarity bound and certainly lie far outside the 
applicability of ordinary leading-twist pQCD. Thus, we expect that for real photons with 
astronomically large energies, a sizable contribution to the total 'yN cross section may arise 
from small size dipoles whose amplitudes in impact parameter space are close to the the 
saturation limit for a wide range of impact parameters. To address quantitatively the issue 
of the real photon - nucleon or the real photon - light nucleus interaction at cosmic ray 
energies we need to use a model which naturally matches with pQCD for small dipole sizes 
as it is the small size dipoles which will play the crucial role in the analysis of maximal 
possible growth of the cross section. Hence we will use the MFGS model. The consistency 
of a particular dipole model can be checked by considering both the real photon limit and 
the limit of extremely high photon energies where the total cross section should start to 
exhibit behavior consistent with S-channel unitarity. We will demonstrate this consistency 
within the MFGS model in the present paper. 

The paper is organized as follows: I 11 Sec. [TT] we discuss the modifications of the MFGS 
model necessary to apply it to real photon - nucleon scattering. In Sec. II I II we evaluate the 
rate of growth of the 7 -proton cross section at UHE and we indicate the advantages and 
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limitations of the current approach. We estimate, within the model, the fraction of the total 
cross section due to the diffractive processes and find it to be quite close to the black disk 
limit of 50%. We also demonstrate that nearly 25% of the total cross section at the highest 
energies comes from production of leading charm (the interaction of the photon in the cc 
component at photon energies of 10 12 GeV). For a target 12 C nucleus, this ratio rises to 
around 30% due to differences in the degree of shadowing for charm and light quark dipoles. 
In Sec. [TV] we extend the results of Sec. nn to the case of a 12 C target. We find that blindly 
extending the procedure we used for the proton target (imposing an upper limit of unity 
on the nuclear profile function) leads to the paradoxical result that the ratio of the nuclear 
and nucleon cross sections exceed the total number of nucleons, A. We explain that this is 
a consequence of the small thickness of the target. The main result of our investigations 
is the general observation that taming is necessary for the elementary cross sections of the 
hadronic subprocesses in y-N scattering, and that, due to the large amount of diffractive 
scattering for the hadronic subprocesses in y-N scattering, there is a significant amount of 
shadowing in y-A scattering. We develop and we summarize our results in the conclusion. 
In App. El we indicate the parameterization that we used to extrapolate the gluon parton 
density to extremely small x, and in App. El we give a general description of the energy 
dependence of UHE cross sections based on the properties of the photon wave function and 
energy dependence of the dipole - nucleon interaction. 


II. THE PHOTO-PRODUCTION LIMIT 


In the perturbative QCD dipole model, the total (virtual) photon - nucleon cross section 
due to the interaction of the small size configurations with a target is written as the convo¬ 
lution product of a basic perturbative cross section for the interaction of the hadronic Fock 
component of the photon with the transverse and longitudinal light-cone wave functions of 
the photon, -0 LiT (,y d), see e.g. Refs. ( ul I t3]): 


a, 


n T N L (x,Q 2 )= dz d 2 d \i) T , L (z,d)\ 2 a tot (d,x') 


( 2 ) 


If the 4-momentum of the incident photon is q , then, in the conventional notation, Q 2 = — q 2 . 
Of course, the longitudinal component of the photon wave function vanishes in the limit that 
the photon is real, Q 2 —> 0, and the transverse component becomes a unique function of the 
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quark momentum fraction, z, and the hadronic size, d. The transversely polarized photon’s 
light-cone wave function is, 


Uf 


\^ T (z,d)\ 2 = -^a e . m . J^eJ[(z 2 + (1 - z) 2 )e 2 A' 1 2 (ed)+ mj/l 0 2 (ed)], 


(3) 


9=1 


where the sum is over rif active quark flavors, rn q is the quark mass, the K's are the 
modified Bessel functions of the second kind, e q is the fractional charge of quark q, and 


e 2 = Q~(z( 1 — z)) + m 2 . (See Ref. j22j for discussion of appropriately dealing with quark 
masses.) For UHE photons, Eq. © will include a term for the light quarks (rn q ss .3 GeV) 
as well as a term for charm (rri q « 1.5 GeV). More massive quarks are strongly suppressed 
by the light-cone wave function and are neglected in the present analysis (but see Sec. IIIII 
for more discussion of the heavier quarks). The energy dependence in Eq. © enters through 
the Bjorken-x variable commonly used in deep inelastic scattering: 

_ Q 2 

2 P-q 

where P is the 4-momentum of the target nucleon. The MFGS model specifies the small 
size and large size cross sections and a scheme for interpolating between the two. In the 


original formulation of the MFGS model 


22 ], x and Q 2 are taken as input for the calculation 
I that the value of Bjorken-x, x', used in the 


of a particular cross section. Note in Eq. 
basic dip ole-nucleon cross section is not the same as the external value of Bjorken-x. A 
relationship between x and x 1 must be supplied by the specific model, and for the MFGS 
model x' is an effective average Bjorken-x for the dip ole-nucleon scattering subprocess, 

, Q 2 4m;;\ 0.75A 


x = 


2 P-q 


1 | 4m 9 

Q 2 


1 + 


(4) 


d 2 {Q 2 + 4m^)_ 

We will now give a brief review of where this expression comes from. The method of 
calculating x' used in Eq. © is derived for the small size components of the photon where 
the photon fluctuates into a qq pair which then exchanges a single gluon with the target. It 
can be verified that the important values of x 1 are centered around x' ~ 1.75x in the usual 
deep inelastic Bjorken scaling regime. For a detailed derivation of the relationship between 


x and x', see Ref. 


By ignoring aligned jet configurations (z —> 1,0), we can write the 


expression for the invariant mass of the produced qq system as, 


M qq ~ 4(mJ + ki). 


( 5 ) 









Using Eq. © along inside the exact expression for x', 


x = 


Q 2 + M l 

s + Q 2 


(6) 


we arrive at the expression, 




1 + 


ki 


(Q 2 + 4 m 2 ) 


(7) 


The typical numerical value kj_ contributing to the integral should be similar to the average 
(ki) ~ X/d 2 that we use in sampling the gluon distribution. In order to reproduce the 
x’ ~ 1.75a; behavior in the Bjorken limit, we use ky = .75A /d 2 . This reproduces Eq. Q 
The model is well-defined (finite) in the photo-production limit, Q 2 —> 0; but the external 
value of x vanishes for all energies when Q 2 = 0 and is clearly not appropriate as input to 
the basic cross section. Rather, we would like to fix Q 2 = 0 and specify £7,. To this end, we 

0 : 


note how the original value of x' used in Ref. 112] 

4 d 2 m 2 ) 


, Q 2 ^0 

X = 


2 _ n 4m 2 


behaves for Q 2 
0.75A \ 


( 8 ) 


2 P-q 

The constant, A, relates the size of the hadronic configuration to the virtuality, Q 2 , of a 
particular quantum fluctuation of the photon through the relation, Q 2 = X/d 2 . Note that, 
even when x —> 0, the effective value, x', is large when d —> 0. We emphasize that, although 
the appearance of Eq. (JHJ) seems at hrst glance to be somewhat ad hoc for use at UHE, it 
is sufficient for our present purposes to simply note that x' is directly proportional to x in 
this limit and inversely proportional to d 2 . Estimates of the sensitivity to other factors in 
expression dHJ) can be obtained by varying the quark mass and the parameter A. 

In the original formulation of the MFGS model, the large size configurations were charac¬ 


terized by growth with energy that mimicked the In 


2 ( W 2 ' 

-w%' 


growth of the pion-nucleon total 


cross section with 


assumed in Ref. 0 that for large sizes, 


Vq = 400 GeV 2 (or x 0 ~ 0.01 for DIS kinematics). That is, it was 


cr(x', Q 2 ) = <j nN (x/ Q 2 ) = 23.78 mb ) 

\ X ' 


0.08 


This behavior for the pion cross section was extracted from data in Ref. 


(9) 


28j ]. Furthermore, 




with a' = .25 GeV 2 and cL = .65 fm. 


soft Pomeron exchange leads to a factor of e 
in the scattering amplitude. In the case of the virtual photon discussed in Refs. j22|,|l|, the 
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external value of x was used for the large size behavior. Clearly, this is inappropriate in the 
photo-production limit where x = 0. 

In fact, we will now argue that the appropriate value of x to use for the large size 
configurations is the same x' (Eq. Q) that was used for the small size behavior in Ref. 22]. 
To see this, consider the large W 2 , fixed q 2 , limit of the pion-nucleon scattering cross-section 
(considering, for the moment, q to be the pion 4-momentum): 


W 2 = (q + P) 2 

= M 2 + 2P ■ q + q 2 
*S° 2 P ■ q. 


( 10 ) 


Then, 


In 


W 2 

W 2 


= In 


= In 


2 P-q 
2P-qo 
x 0 


X 


+ In 


91 

Q 2 0 


(ii) 


Where Q 0 is defined so that x 0 = Ql/(2P ■ q 0 ) = .01. If we regard the pion mass as a 
particular value for the photon virtuality, then we see that Eq. m generalizes to any small 
photon virtuality Q 2 and small Bjorken scaling variable, x. On the other hand, in the limit 
of small Q 2 and non-vanishing size, d, Eqs. 0] |H1 and yield, 


In 


^ 1 = hr 
x' 


2 P-q 


= In - 


( 12 ) 

(13) 


where x r 0 corresponds to taking Q 


2P ■ q 0/ 

fW 2 ' 

w. 

Q 0 . In other words, to truly mimic the behavior of the 
pion-nucleon cross section, one should use the effective Bjorken-x that was used in sampling 
the gluon distribution. Comparing Eqs. and E3 we also have 

ln (5) - 111 (?) =111 (|) ■ (14) 

Thus, if one uses the external values of x in Eq. © then one over-estimates the large size 
cross section by a factor of ( Qq/Q 2 )' 0S which diverges in the limit of a real photon. Finally, 
we can determine the value of x' Q by noting that Wq = 400 GeV 2 implies Ql = 4 GeV 2 . 
Therefore, from Eq. (JHJ), we have, 

.75A 


x' 0 = x 0 (l + m 2 q ) ^1 + 


4d 2 (l + m 2 ) 


(15) 
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The second term, m 2 , in parentheses in Ref. E3is implicitly divided by 1 GeV 2 so that it is 
unitless. Notice that the difference between Xo and x' Q is only significant for small d or large 
m q . The conclusion of this section is that it is appropriate to use the effective x' given in 
Eq. (0) for all values of Q 2 and that Eq. (j3J) must be used in the photo-production limit. 
Furthermore, x' is calculated unambiguously from Eq. 0 and E 7 and with the condition 
that the target nucleon is at rest so that P ■ q = ME 1 . For the rest of this paper, we assume 
that the target nucleon is at rest and we specify E 7 as input for the dipole model of the real 
photon. 


III. GROWTH OF THE CROSS-SECTION AT VERY HIGH ENERGIES 


W 


len one considers the profile function as in Ref. 


231 ]. or the S-matrix as is done in 


Ref. 3], one finds that the unitarity constraint is usually violated at very low values of 
Bjorken-a; (Ultra-High Energy). This sort of behavior is evidence of breaking of the leading 
twist approximation and indicates the onset of qualitatively new QCD phenomena. However, 
the small size configurations of the UHE real photon wildly violate the unitarity constraint 
due to the rapid growth of the perturbative expression, 


7r 


VpQCDid,* 1 ) = ^-d 2 a s {Q 2 )x'g N (x', Q 2 ), 


(16) 


at small x' even when the increase of the radius of the interaction with energy is taken into 
account. Therefore, it becomes necessary to introduce some new assumptions about the 
unitarity violating components of the total cross section in order to make some sense of the 
UHE behavior. We choose to adhere to the usual assumption that at UHE, in the region 
where small size perturbative methods begin to break down, the cross section is at (or near) 
the limit set by unitarity. This assumption is supported by our studies of the amplitudes 
of the dipole - nucleon interaction at HERA energies j^. For reference, the hadron-proton 
amplitude that was used in Ref. 123] and which we use in this paper is given by, 

1 1 


with, 


Thv'('S)t) is O’tot 


M 2 (d 2 ) = 


a 1 ^ In (HO 


(1 — t/M 2 (d 2 )) 2 1 — td 2 /d 2 ml 


m t 


m 


(m 


i ~ m l )S , d<d n 

U '7T 

, otherwise 


(17) 


(18) 
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(See Ref. Q] for the detailed procedure for obtaining this.) For these equations, rn\ m 
1.1 GeV 2 , ml ~ 0.7 GeV 2 and m 2 ~ 0.6 GeV 2 . The typical size of the pion is d n & 0.65 fm, 
a' is 0.25 GeV -2 , and xo = .01. We use Eq. (fT7l) to provide continuity with our previous 
publications while noting that the details are not important for our current purposes; any 
reasonable model of the t-dependence that matches correctly to lower energy behavior is 
sufficient. At UHE, the dominant effect is from the diffusion factors. By making a large 
estimate for the amount of diffusion (a '), we over estimate the amount of spreading in impact 
parameter space. Within our approach to unitarity taming, as will become clear, this over 
estimate of the diffusion is consistent with the goal of obtaining an upper limit to the growth 
of the cross section. 

Very large size configurations tend not to violate the unitarity constraint in a gross way - 
the total cross section for large size configurations increases slowly with energy, but there is 
suppression of the impact parameter at small b with increasing energy due to diffusion con- 




does violate 


sistent with the Donnachie-Landshoff soft Pomeron j29| . In fact, the fit of 
S-channel unitarity for LHC energies and above. However, this leads to a very small effect 
for the total cross sections for the energies we discuss here. We can enforce the unitarity 
assumption about the high energy behavior if we re-calculate the total photo-production 
cross section by using the basic cross section of the MFGS model with linear interpolation 
everywhere except where the unitarity limit is violated. When the unitarity limit is violated, 
we set the cross section equal to the maximum value allowed by unitarity. More specifically, 
we first decompose the the basic cross section in terms of the impact parameter represen¬ 
tation of the amplitude (usually called the profile function) using the optical theorem. Let 
l) be the hadron-nucleon amplitude and let l be the 4-momentum exchanged in the 
subprocess. Then, we calculate, 


=21 dz I d‘d \ii T (z,dW I d 1 bT(E 1 ,b,d), 


(19) 


where r(E 7 ,6, d) is the profile function for a configuration of size d. The hadronic profile 
function that we use is given by the usual definition, 


T( E „b,d) = 2 1 J d 2 le ihl F hN (Ey, /), (20) 

so long as it is less than one. If it exceeds unity, then we explicitly reset T(Ey,b, d) = 
1.0. Note that in the limit of d —> 0, both x’ and the photon wavefunction become mass 
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Light Quark Contribution 



d (fm) 

FIG. 1: The distribution of the integrand in Eq. 0 over hadronic sizes. Here, the unitarity 
constraint, T < 1, is explicitly enforced. Nevertheless, the distribution becomes sharply peaked 
around small hadronic sizes for UHE photons. 

independent. As we shall see, the distribution of hadronic sizes involved in the interaction 
becomes more and more sharply peaked around small sizes as Ej —> oo. Thus, in the UHE 
limit, the sum in Eq. © will contain significant contributions from all massive strongly 
interacting particles. The energies that we consider in this paper are still not high enough 
to include all heavy particles, but, as we will discuss, there is a significant charm contribution. 
We include the charm contribution in the calculation of the cross section and note that the 
resulting hadronic showers will likely consist of a significant number of charmed particles like 
D-mesons. The contribution from the bottom quark will be suppressed relative to the charm 
due to its larger mass and smaller electromagnetic coupling. We neglect the contribution 
from bottom and all heavier particles. 

In our calculations, we use CTEQ5L gluon distributions js(|. We note that, since sig¬ 
nificant sections of the matching region will violate the unitarity constraint, then trying to 
achieve an extremely smooth interpolation is an arbitrary modification to the model which 
achieves no genuine improvement. For simplicity, therefore, we use the linear interpolation 
in the MFGS model. The plots in Figs. |T| and [21 demonstrate the resulting distribution of 
hadronic sizes in the photon. Naturally, the peak at small sizes becomes sharper in the UHE 
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Charm Quark Contribution 



d (fm) 

FIG. 2: This is the distribution of the charm contribution normalized to one. Comparing with 
Fig. ^ we see the suppression of massive quark contributions. Comparing this with Fig. ^ we 
see that the distribution for heavier quarks is more sharply peaked around small transverse sizes 
relative to lighter quarks. 


limit. 


Before leaving the subject of t-dependence in the UHE photon, we recall that in the 


original MFGS model of the t-dependence, discussed in Ref. 


23], the diffusion of the small 


2c 


size qq pairs was neglected. However, in a recent overview 1311 of the behavior of hadronic 
cross sections at UHE, the maximum diffusion of small size configurations at UHE was 
considered. In order to correct for the small size diffusion, we slightly modify the diffusion 
factor in the amplitude of Ref. to the form, 

F F (t,x) = e“'U)f ln Uo/x), ( 21 ) 

where a'(d ) = .25(1—0.5e -18cfl ) GeV~ 2 . For sizes, d, greater than the pion size (d n ~ .65 fm), 
the value of a' quickly approaches the usual slopefcd = .25 GeV~ 2 ) in accordance with the 
Donnachie-Landshoff soft Pomeron as in Ref. |23(. However, unlike the Regge slope in 
Ref. Q which vanishes at small sizes, the value of a' approaches .125 GeV -2 for d < .2 fm, 
which is the diffusion rate for small size dipoles determined recently in Ref. laiJ . The use 
of the exponential function ensures that the interpolation through the transition region is 
continuous, and that a' approaches .25 GeV -2 rapidly for d > .65 fm and .125 GeV -2 for 
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d = 0.1 (fm) 


d = 0.3 (fm) 


d = 0.6 (fm) 


1 

0.8 

0.6 

L, ^ 

0.4 

0.2 

0 

0 1 2 0 1 2 0 1 2 
b (fm) b (fm) b (fm) 

FIG. 3: Samples of the profile function for the hN interaction for real photon energies of 10 4 GeV 
(dashed line), 10 6 GeV (dotted line), 10 8 (dot-dashed), and 10 11 GeV (solid line) for a range of 
hadronic sizes. 

d < .2 fm. Note that at very high energies we should take into account that the Fourier 
transform of F iP (i, x) should contain a tail oc exp — fib for some mass scale, /i. However, in 
the energy range discussed in this paper it is a small effect and hence we neglect it. 

Samples of the profile function obtained when we use the above f-dependence and the 
above procedure for taming the unitarity violations are shown in Fig. 01 Note the extremely 
rapid growth with energy at sizes of d — .1 fm. 
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y-proton cross section 




FIG. 4: Growth of the total photon nucleon cross section for the range of energies from E~, = 
10 3 GeV to = 10 12 GeV. On the x-axis of the bottom panel, ln 3 (^-) (Eq = 1.0 GeV) is plotted 
to allow for easy comparison with the expected energy dependence at UHE (see App. 0. The 
lowest three curves in each panel are the result of using the MFGS model with variations in model 
dependent parameters to demonstrate numerical sensitivity (see text). For comparison, the upper 
curve shows the two Pomeron model of Donnachie and Landshoff Ref. Q and the dotted curve 
shows the model of Shoshi et. al.Ref. 



The procedure described above may be regarded as placing an upper limit on the growth 
of the cross section since we have taken the maximal contribution that does not violate 
unitarity. The resulting nucleon cross section has been plotted in Fig. 0 To test the 
numerical sensitivity to a variation in the upper limit of the profile function, we have included 
the result of placing the upper limit of the profile function at 0.8 rather than at 1.0. In 
the studies jljj, | 22 j the matching parameter, A was estimated based on the analysis of 
the expressions for <jl(x,Q 2 ) to be of order 10. A later analysis of the J/-0 production Q 
suggested that a better description of the cross section for the intermediate 0.5 > d > 0.3/m 
is given by A ~ 4 while the cross section in the perturbative region depends very weakly on 
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A. We therefore include in Fig. H] the result of using A = 10 to test the sensitivity to the 
matching ansatz. It is evident that variations in these parameters make a small reduction 
in the upper limit. For the rest of this paper, we assume that A = 4 and that the unitarity 
constraint implies, < 1 . 

For comparison, the Donnachie-Landshoff two-Pomeron model 
constraint is also shown in Fig. 01 Yet another approach was used for 
The result of this calculation is shown as the dotted curve in Fig. 0] In the model of 
Ref. ( 25 )], the large size components are modeled by non-perturbative QCD techniques, but 
the small size components use a Donnachie-Landshoff two Pomeron model with parameters 
somewhat different from the original model. While the method used in Ref. |25| imposes 
impact parameter space unitarity in the dipole picture, it does not account for the non- 
normalizability of the photon wavefunction. Instead, a unitarity constraint is imposed upon 
the S'-matrix for 7 -proton scattering, and is thus less restrictive than our approach which 
applies unitarity constraints to the profile function for individual hadronic configurations. 
Our approach thus restricts the growth of the cross section more than either of the two 
above scenarios as can be seen in Fig. 0 

It has been observed that agreement with data is improved if one uses MS-bar NLO pdfs 
in very low-x DIS experiments. Our main concern is that the rate of increase described by 
the interpolation that we obtained in Eos. IA1IIA3I of the appendix is not drastically modified 
by the inclusion of higher order effects. In particular, the rate of growth, and the order 
of magnitude of higher order corrections should not be drastically altered from what one 
obtains at leading order. To argue that this is the case, we notice that the LO and NLO 
parton densities grow at nearly the same rate at very high energies. This is demonstrated 
for a typical small configuration size in Fig. 0 where we compare the LO and NLO CTEQ5 
parton densities lowest values of x where parameterizations exist. The main difference 
between the two parton densities is that, in the high energy limit, the leading order gluon 
density is nearly a constant factor larger than the NLO pdf. (This makes sense when one 
considers the color factor 9/4 difference between quark-antiquark dipoles, and gluon dipoles 
and the need to include qqg dipoles in a consistent NLO formulation of the dipole model.) 
Indeed, above about = 20000 GeV in the plot of Fig. O the LO curve is shifted by about 
a constant factor of 9/4 upward from the NLO curve, and the energy dependence used (see 
Eq. O of App. 0 ) describes both curves with an accuracy of < 10 % for the energy range 


1 32] with no unitarit 


alculations in Ref. 


taritv 
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FIG. 5: Plot comparing the the growth of LO and NLO pdfs at high energy for the small dipole 
size, d = .05 fm. 


Ey > 20, 000 GeV. In the spirit of obtaining an upper limit on the growth of the cross 
section, we continue to used the LO pQCD result. Since we checked that our LO inspired 
parameterization of a(d,x') describes the data well down to x ~ 10 -3 the difference in the 
normalization is likely to be partially absorbed in the definition of the cross section. Note 
also that the recent studies of the small x behavior of the the gluon densities indicate that 
the NLO approximation is close to the resumed result down to x ~ ICO', see a review in 


EL 


Before leaving the discussion of the nucleon target, we note several advantages and disad¬ 
vantages of the current approach: First, we note that Fig. Elis consistent with the rate growth 
of the cross section for a real photon at UHE which we find in App. El to be r-u In \Ey/E 0 ), 
where E 0 is of order 1 GeV (see e.g. f29j for a review of the Froissart bound for the hadron 
- hadron scattering cross section, See App. |B] for a generalized discussion of energy depen¬ 
dence for real photons.) We are also able to analyze the contribution from different flavors 
separately. In Fig. El we see that the fraction of the total cross section due to charmed 
particles rises to nearly 25% at 10 12 GeV [^7|. This indicates that a large fraction of showers 
initiated by UHE photons should contain a pair of leading charmed particles. It is worth 
investigating whether such showers have a substantially longer penetration depth in the 
atmosphere and could be separated by the Auger detector. 

Furthermore, by recalling the relationship between the diffractive component of the basic 
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fractional contribution of charm 



FIG. 6 : The fraction of the total hadronic 7 -proton cross section due to charm qq pairs. 

diffractive fraction of cross section 



FIG. 7: The fraction of the total hadronic 7 -proton cross section due to diffractive reactions of 
the hadronic component. The large fraction of the cross section due to diffractive scattering at 
Ery = 10 12 GeV indicates relative nearness to the black disk limit. 

cross section, &diff and the hadronic profile function, 

CTMff = f d 2 b\r h (b)\ 2 , (22) 

we may also separate the fraction of the total hadronic profile function due to diffractive 
scattering. This is shown in Fig. [71 which gives a integrated measure of the proximity to the 
black disk limit. In the black disk limit, diffractive scattering accounts for exactly half of the 
total cross section. Therefore, the fact that, as Fig. [3 shows, the fraction of the cross section 
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due to diffractive scattering is around .35 at UHE indicates that diffractive scattering plays 
a significant role and that there will be large shadowing in nuclei. 

Since the cross section grows extremely quickly at UHE, then the unitarity limit is sat¬ 
urated even at values of dipole size around d — 0.1 fm. Figure [I] demonstrates that the 
largest contributions to the total photo-production cross section come from regions around 
d « .1 fm and from the transition region. Sizes smaller than this contribute very little 
to the total cross section. Thus, pQCD provide very little detailed information, since all 
models of the basic cross section (at small sizes) which violate the unitarity constraint at 
d « 0.1 fm will give very similar results. Furthermore, since the gluon distribution rises 
very sharply between d — 0 fm and d — 0.1 fm at very high energies, then the calculation 
in Eq. m becomes more sensitive to how the gluon distribution is sampled. Thus, there is 
more sensitivity to the parameter, A, used to relate the hadronic size to the hardness of the 
interaction. As seen in Fig. 01 if A = 10, then the cross section is suppressed by around ten 
percent from its value when A = 4. 

The largest source of theoretical uncertainty comes from the region of large hadronic 
sizes. Though the behavior of large configurations can be reasonably expected to follow 
pion behavior at accelerator energies, we do not have any experimental data for > 10 6 GeV 


hadrons with which to model these extremely high energy Fock states. Moreover, we have 
so far been associating each hadronic Fock state with a particular size. It may be that as the 
energy of the photon increases, a large number of hadronic Fock states (perhaps multiple 
pion states) may be associated with a single size. Moreover, contributions from large impact 
parameters become significant for extremely energetic photons. Thus, predictions become 
sensitive to how the model handles the f-dependence of the amplitude. The current model 
is based on the assumption that the typical t-dependence for low energies continues into 


the UHE regime. Note, however, that in the case of the pp scattering the analysis of [311 
indicates that, though the black disk limit leads to the slope B oc In 2 s, this is a very small 
correction for the energies discussed here. Finally, we want to emphasize that we allow the 
impact parameter amplitude to approach T = 1 without a slowdown at T ~ 0.5 as happens 
in many other dipole models where eikonal type parameterizations of the dipole-nucleon 
cross section are used (see, for example, [35]). One again it is consistent with our aim of 
yielding a conservative upper bound. 
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IV. NUCLEAR TARGETS 


Since we are interested in the interaction of UHE photons with atmospheric nuclear 
targets, we now go on to investigate the growth of the UHE cross section for the case 
of a real photon scattering from a nuclear target. The steps apply to any of light nuclei 
constituting the atmosphere, but we use 12 C for the purpose of demonstration since the 12 C 
nucleus has the approximate number of nucleons for a typical atmospheric nucleus. 

At first glance it seems natural to repeat the procedure we followed for the proton case in 
Sec. lIIIl with the profile function given by a nuclear shape, and cross section for the interaction 
of the small dipoles as given by Eq. (flU . However we found that if we follow this procedure 
we end up with the obviously wrong result that in the UHE limit, this approach quickly leads 
to the situation that ct 7 a >> which is physically unreasonable. The reason for this 

becomes clear if we visualize the relationship between the proton PDF and the nuclear PDF. 
The unitarity condition is not sensitive to effects of transverse correlations of the partons. 
The unitarity constraint would tame the dynamics if one could assume that the nucleus is a 
perfectly homogeneous distribution of nuclear matter. At high energies, the disk of nuclear 
matter “seen” by the incident hadronic configuration blackens as represented schematically 
in Fig. [HI Any inhomogeneity in the distribution of nucleons is accounted for at low energies 
in the grayness of the nuclear disk without yielding a quantitative difference in the total 
cross section. However, the actual distribution of nuclear matter in light nuclei “seen” 
by the incident dipole probably looks more like that of Fig. |H1- a collection concentrated 
regions of nuclear matter which individually grow black in the high energy limit but far apart 
from each other transversely. If the nuclear system is dilute and nucleons do not overlap 
transversely, the use of the T(b) < 1 condition becomes insufficient. Thus, we will certainly 
find that a 1 a >> Aa^ if we assume that both the disk of the individual nucleon and the 
disk of the nucleus grow black in the high energy limit. That is, the cross section resulting 
from Fig. Mb) is certainly larger than the sum of the cross sections from each of the blackened 
nucleons seen in Fig.lHJb). The simplest illustrative example would be to consider scattering 
off the deuteron - in this case neglect of the cluster structure of the system would grossly 
overestimate the maximal cross section for the interaction of this system with a small dipole. 
It is worth emphasizing that all these considerations are valid for the light nuclei and are 
likely to be a correction for the scattering off sufficiently heavy nuclei. A more meaningful 
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FIG. 8: The nuclear disk as “seen” by the incident hadronic configuration in the simple homoge¬ 
neous model of the unintegrated nuclear PDFs. The level of absorption by the disk is indicated 
by the level of grayness. At low energies, (a), the disk is weakly absorbing and homogeneous. At 
very high energies in (b) the disk becomes black and is thus totally absorbing. 



FIG. 9: A more accurate way to visualize what the incident hadronic configuration probably 
actually “sees”. The nucleus consists nucleons separated over a large distance with concentrated 
nuclear matter in (a). At very high energies (b), each of the nucleons becomes totally absorbing 
(black). 


upper estimate of the cross section of scattering off light nuclei is, therefore, to take into 


account first the taming of the elementary cross sec 


theory of nuclear shadowing due to diffraction 


36 


ions and next the Glauber - Gribov 


37] which does not rely on the twist 


decomposition of the cross section. Since we observed that diffraction constitutes a large 
fraction of the total cross section we expect that a large shadowing effect will emerge with 
growing energy. Consequently, our result will automatically be consistent with S-channel 
unitarity. Hence, we take the usual approach to nuclear scattering when the product of the 
nuclear optical density with the cross section is small: Aa^jyT(b) < 1. Following in the spirit 
of the_treatment of hadronic fluctuations in the Good and Walker picture as presented in 
e.g. 38], we consider the states of the incident photon to be a linear combination of nearly 
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“frozen” hadronic states that do not mix with one another during their passage through the 
target. Each of these states scatters from the target with a cross section (&h) = dh- Let us 
write the total photon-nucleus cross section as, 

a lA = Aa h - A a. (23) 


A standard result of the Glauber-Gribov theory in the language of hadronic fluctuations is 
that the full shadowing correction, Act, can be written as, 

A a=j J d 2 hT 2 {b)(a 2 h )e-^ A{dh)T{ - b) . (24) 

Where T(b) is the nuclear optical density normalized to unity. Equation (1241) is an approx¬ 
imate formula valid for the case of small fluctuations or small nuclear thickness (the later 
is true in our case). We will find a convenient expression for A e ff/A if we expand this 
expression using the argument of the exponent as a small parameter: 

a lA = Aa h - J d 2 b T 2 (b ) + • • • (25) 

If A e ff is defined by the relation, a yA = A e fftj lN , then dividing Eq. ( 1251 ) by Aa' yN gives, 

; W = i -LL/d 2b r 2( t ) + ".. (26) 


Define an effective cross section, 

a eff = 

where by definition, 


(d 2 ) 167T dd^ if 


a 


7 AT al N 


t =0 


(a 2 ) = I dz I d 2 d\'i/j T (z : d)\ 2 dl{d,x'). 


Then we can write, 


(27) 


(28) 


Aeff <Teff - f d 2 bT 2 (b) + ■ ■ ■ 

A aeff 

I d 2 b ( T(b)a e „ - % T 2 (b) + 

a eff 

Considering the first two terms in Eq. (j2HI) as the first terms in a power series expansion 
we may write, 


(29) 


m, 


A eff _2f d 2 b (1 - e~^ff T(b) ) 


A 


a eff 


(30) 
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Effective Profile Functions 


Carbon-12 



FIG. 10: Sample of the effective profile functions (see Eq. (EH) for a real photon on a 12 C target. 

If we identify the effective profile function as, 

T eff (b) = l- e -^ T(b) , (31) 

then we see that both the effective profile function and the shadowing ratio are less than 
unity by construction. Furthermore, because A e ff must be less than A (for scattering in any 
range of impact parameters), then we may regard F e ff(b) < 1 as our unitarity condition for 
the nucleus. 

(Note that the diffractive components have mass squared proportional to 1/d 2 and there¬ 
fore correspond to rather large hadron multiplicities.) The result of evaluating Eq. © 
and solving for the 7 A cross section is shown in Fig. One can see that the energy 

dependence of the y- 12 C scattering cross section is substantially weaker than for the proton 
target due to nuclear shadowing, though the increase of the cross section as compared to 
the energies studied experimentally is still large. Since there is less shadowing for the case 
of incident charm dipoles, then our analysis indicates that there is a larger fraction of the 
total photon- 12 C cross section that arises due to charm dipoles than in the case of a proton 
target. This is shown in Fig. EH In the case of a target 12 C nucleus, the fraction of the 
total cross section due to charm dipoles is around 30%. The large amount of shadowing 
that we find already has implications for ener gies around ~ 100 TeV which are relevant to 
a number of current cosmic ray experiments (22|. In addition, the forthcoming studies of 
ultra-peripheral heavy ion collisions at the LHC would allow, to some extent, a check of 
our predictions by measuring shadowing for photon-heavy ion interactions for the range of 
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y + 12 C cross section 



Shadowing Ratio 



FIG. 11: The upper panel shows the dependence of the 7 - 12 C cross section on the incident photon 
energy. The lower panel shows the dependence of the shadowing ration ct^a/{Act^n)- 


values for from 1000 GeV to 2000 GeV. Within our model, we find that A e ff/A « .3 for 
a/s = 100 GeV and A e ff/A ps .2 for a/s = 2000 GeV with A = 220. In order to allow for a 
simple extrapolation from the shadowing ratio, ( A e ff/A)c , for Carbon to other nuclei with 
masses typical of atmospheric atoms we use the function, 


A 


eff 


A 


A 


eff 


A 


C 


12 


(32) 


where we then determine n for a a set of fixed photon energies and for the range of atomic 
masses, 12 < A < 16. As a sample, we list the following: For E 1 = 10 12 GeV, we find 


n = —.41; for E~ = 10 9 GeV, we find n = —.35; and for E 1 = 10 6 GeV, we find n = —.3 


— 1 n6 


V. CONCLUSION 

In this paper we have investigated the photon-proton/nucleus cross section in the range 
of energies from 10 3 GeV to 10 12 GeV. Figure 0 demonstrates that the total cross section 
rises by about a factor of 12 , but that there is a significant amount of uncertainty involved 
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Fraction of y-C 12 Due to Charm Dipoles 



FIG. 12: The fraction of the total 7 - 12 C cross section due to charm dipoles. 


when the details of the model are varied. This result, however, gives us a very reasonable 
upper limit on the cross section since we have consistently taken the maximum cross section 
allowed by unitarity. The cross section varies approximately linearly with ln 3 (^-) (see 
App. El). (We use E 0 = 1.0 GeV which is consistent with hadronic sizes; see, for example, 
page 18 of Ref. jstj])- 

The solid curve in Fig. 0 provides a reasonable estimate (or if we prefer to take a more 
cautious attitude, an upper limit) to the 7 IV cross-section at extremely high photon energies. 
We note that the dipole approach is consistent with the direct extrapolation of the photo- 
nuclear cross section |40j and with a model based on unitarity in the f-channel 411, but 


not with the model of |32{] which uses a two-Pomeron approach and places no unitarity 
constraint on the growth of the cross section. There is also disagreement with the model of 


• 0 


Ref. 123 which uses the Donnachic-Landshoff two-Pomeron approach, but applies a unitarity 
constraint to the 5-matrix for 7 -proton scattering rather than to the profile function of 
individual hadronic components. 

We have shown in Fig. El that as much as 25% of the cross section may be due to charmed 
mesons for the case of a target proton, and this fraction rises to 30% when we consider 
a target 12 C nucleus. For the 12 C target, we use a direct application of the usual Gribov- 
Glauber theory, and we find the there is a large amount of shadowing increasing with energy. 
It is worth noting here that though our results for the elementary yp cross section are rather 


close to the results of 


for \fs < 10 3 GeV considered in this paper using a generalized 
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vector c 
of Ref. 


ominance model with a point-like component in the photon wave function, the model 


42] leads to a nuclear shadowing effect which is practically energy independent. 


This is consistent with the dipole-proton interaction approach to the black disk limit. The 
resulting cross section is shown in Fig.^Jwhich indicates a rise in the cross section of about 
a factor of 7 when the energy increases from 10 3 GeV to 10 12 GeV. The relevance of these 
observations is in the characterization of atmospheric showers induced by UHE neutrinos 
and super-GZK cosmic rays, where upper limits on the allowed growth of the photon-nucleus 
cross section are needed. 
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APPENDIX A: INTERPOLATION TO VERY SMALL BJORKEN-X 

In order to further extrapolate the basic, small size, cross section to extremely small values 
of x, we make a fit to the CTEQ5L gluon distribution in the region of lowest x (10 -5 > x > 
ICE 4 ) where the parameterizations exist. We find that the following interpolation agrees to 
within a few percent over the range of small-x and for .01 fm > d > .2 fin: 


xgjy(x,d) = a(d)x c ^ 

(Al) 

c(d) = —.28 d~' u 

(A2) 

a(d) = 3.9 - 13.9d + 20.1d 2 + .5 In d. 

(A3) 


27 






APPENDIX B: LIMITS ON ENERGY DEPENDENCE AT ULTRA-HIGH EN¬ 


ERGIES. 


It is possible to obtain the general energy dependence, a lh ~ In 3 for the 7 -hadron 
cross section in the UHE limit within the dipole model. Here we give a general proof based 
only on the following assumptions about the UHE limit: 

• The dipole model for a finite number of active quark flavors holds for the real photon 
in the UHE limit in the sense that at a fixed energy, the dipole cross section increases 
with d no faster than d 2 . 


• A given finite size hadronic Fock component of the real photon scatters with exactly 
the maximum possible cross section allowed by the unitarity constraint when E 1 —> 00 , 
and the rate of increase of the cross section for each individual Fock component of size, 
d, is limited by the rate of growth in the Froissart limit. 


• The very small size hadronic Fock components of the real photon scatter with a cross 
section whose rate of growth is no faster than a power of x' (or Efj . 

The second bullet above requires some clarification. Usually, the Froissart bound is only 
applied to the interaction between two hadrons rather than to the interaction between a 
wave-packet and a hadron. However, since the derivation of the Froissart bound is based 
on analyticity in the t-channel which leads to the requirement that the amplitude in impact 
parameter space falls of at least as fast as e~ 2m,rb , 2 , 9 ] then the argument works in our case 
as well (49f. 


We are only interested in the variation of the cross section with energy. Therefore we 

will leave out over-all factors in order to simplify the argument. Note that In A can always 

x 0 

be separated into a sum of ln(i? 7 ) and terms that only depend on d. For the rest of this 
section, we will always write In as liUFL) since it is the leading powers of photon energy 
that will interest us. To be concise, the symbol ~ will indicate how a cross section varies 
with photon energy, U 7 , whereas oc will indicate how a cross section varies with hadronic 
size, d. The first of the above assumptions allows us to state that, 


u(U 7 


E 

flavor 


ddd\^{d)\ d(d,x r ) 


(Bl) 



The integral over momentum fraction from Eq. PI) is assumed to be implicit. Also, for the 
rest of this section, the sum over flavors in Eq. ED will be understood and left out. Since 
the energy dependence of the integrand in Eq. m can be understood in the extreme limits 
of d —> 0 or d —> oo, but is model dependent in the intermediate range of d, then let us 
separate Eq. m into the sum of three terms: 


regionl regionl region 3 

'Te * ' ''"7a ^ s 

a(E 1 ) ~ / ddd \i/{d)\ 2 a(d, x') + / ddd \ip{d)\~ a(d, x') + / ddd \i/i(d) | 2 a(d, x') . 

Jo Je J A 

(B2) 

Call the terms in Eq. (El regions 1,2, and 3 respectively. For any given range of photon 
energies, one can choose sufficiently large A and sufficiently small e, that regions 1 and 3 
must give a negligible contribution to the over-all cross section. We will justify this statement 
now. 

First, in region 1, the cross section for the subprocess has the following behavior due to 
the first and last bulleted assumption above: 


d\ (d, x ') ~ d 2 (E 7 )“ , 


(B3) 


where a is some positive real number. Furthermore, the light cone wavefunction of the 
photon gets its energy dependence from the leading behavior of the modified Bessel functions 
in Eq. flu. In the limit of d « l/m q , |?/>(d )| 2 oc 1/d 2 . Hence, in the limit defined by region 
1 , we have the following general energy dependence: 

region l ~ f ddd^d 2 (E 1 ) a ~ e 2 (E 7 ) Q+1 . (B4) 

Jo d 

Next we consider the other extreme: d — > oo. Away from d = 0, Eq. ESI) shows that 
the cross section, d(d, x ') for the subprocess rises very quickly to values that violate the 
unitarity constraint since we are considering UHE photons. At a certain value of d , the 
growth of d(d,x') with d must level out. Within the dipole model, the growth of d-(d,x') is 
flat with respect to variations in the transverse size of the hadronic component in the limit 
that d is large. Call the upper limit of the basic cross section, d max . Also, for d —> oo, the 
Bessel functions give ^(d/f oc e 2r / >q 'J The energy dependence of the large size cross section 
can grow no faster than In 2 E 7 due to the Froissart bound (the second bulleted assumption 
above). Thus, for region 3 we have, 

roo —2 m q d 

region?) ~ / ddd ——— d max ln 2 (E 7 ) ~ e - 2 m 9 A i n 2 E 1 . (B5) 

Ja d 
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For a particular range of photon energies, we may always choose e small enough, and A 
large enough that regions 1 and 2 give a negligible contribution to the total over-all integral 
in Eq. m- From now on, assume that e and A are always chosen, in each energy range, so 
that regions 1 and regions 2 are defined to be negligibly small. Due to the general properties 
of the dipole model there is always a very small contribution from very large hadronic sizes 
( d > A) that grows slowly with energy (~ In 2 Efi), and there is always a small contribution 
from very small hadronic sizes whose contribution may grow very quickly (as a power of 
E y , according to the third bulleted assumption) due to the fact that there will always be 
a contribution from extremely small sizes whose value of cr(d,x') has not yet reached the 
unitarity limit at a given photon energy. We will now consider the rate of growth of the cross 
section that results from assuming that the cross section attains the maximum value allowed 
by unitarity for the largest range of sizes possible within the general constraints of the dipole 
model. Regions 1 and 3 give negligible contributions to the total integral, as discussed above, 
and the values of cr(d, x') will be assumed to saturate the unitarity constraint for all values of 
d outside of range of region 1. This means that for region 2, the basic cross section (denoted 
by a subscript 2 ) has reached the maximum allowed value, cr max , in terms of its growth with 
d, and the rate of growth with E 1 is the maximum allowed by the Froissart bound. The 
cross section appearing in the integrand of region 2 then becomes, 

cr 2 (d, x') ~ cr max In 2 E 1 . (B 6 ) 

However, the requirement that region 2 contains all of the unitarity saturating contribution, 
and that region 1 contains a negligible contribution demands that we allow e to have some 
energy dependence. This is because, as Eq. (ED> shows, the basic cross section at small 
sizes may have as much as a quadratic d-dependence and potentially very rapid energy 
dependence. Therefore, at a small but fixed value of d , the basic cross section quickly rises 
from small values to unitarity violating values with increasing energy. However, region 1 
by definition contains only the suppressed part of the integrand near d = 0, whereas the 
unitarity saturating region should be associated entirely with region 2 . As the energy of the 
photon is increased, therefore, we must continuously redefine region 1 so that the integrand 
of region 1 is confined to a smaller and smaller region around d = 0. This sort of behavior 
does not exist at d > A, because at such large values of hadronic size, the cross section 
only increases as In 2 E 7 and there is almost no variation with hadronic size. Therefore, A 
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is defined without any energy dependence (A may be given weak energy dependence, but 
that will only result in subleading powers of In E 1 in the final result.) Equation (111Ml) tells 
us that the maximum rate at which e may decrease at small d is, 

e(E 7 ) ~£ 7 - (1+q) . (B7) 


We thus write region 2 as, 


region 2 ~ 



cl del ^ @ max In E~, 


~ ln 3 4 (E 7 ) 


(B8) 


Here we have continued to use 1/d 1 2 behavior for the squared photon wavefunction because 
this yields the fastest possible rate of divergence of the wavefunction at the lower end of 
the integral and thus yields the most conservative upper limit on the rate of growth. Notice 
that after having given e energy dependence, Eq. (HU) becomes, 

r e i E ~t) i 

region 1~ / ddd—d 2 \v/ 2 ( K E 1 )^E~^ 1+a \ (B9) 

Jo d 

Thus we have established the behavior of each of the regions in Eq. CD. Regions 1 yields a 
vanishing contribution to the total integral as E 7 —> oo and region 3 has energy dependence 
< In 2 E 7 whereas region 2 has energy dependence, < ln 3 (E 7 ). Taking the leading behavior 
in Eq. CD therefore gives, 


<7(7 hadron —> hadrons) < Constant x In 3 E 7 . 


(BIO) 


Equation (IB 101) applies for each active flavor individually and thus for the sum of flavors. A 
possible way that the rate of growth at ultra-high energies violates Eq. IB 101) in spite of the 
unitarity limit being saturated for each flavor would be for there to be a large proliferation 
of new active flavors at ultra high energies. 
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